Current Research in Structural Biology

Ubiquitin-like protein 3 (UBL3) is a post-translational modifier that sorts proteins into small extracellular vesicles and regulates the trafficking of disease-associated proteins such as -synuclein. The structural and dynamic features of the UBL domain that underlie these functions, however, remain poorly understood. Here we performed in silico structural dynamics analysis of the UBL3 UBL domain using an NMR structure ensemble combined with anisotropic network modeling (ANM) and perturbation response scanning (PRS). Principal component analysis and residue-wise fluctuation analysis consistently revealed high flexibility in the C-terminal region of UBL3. Comparative ANM analysis across 20 ubiquitin-like proteins (UBLs) further showed that C-terminal flexibility is a conserved yet variable property within the UBL family. PRS analysis demonstrated that residues forming the central -helix of the {beta}-grasp fold exert greater dynamic control over collective motions than {beta}-sheet residues. Notably, UBL3 displayed the highest helix/sheet PRS effectiveness ratio among all UBLs analyzed, highlighting the prominent dynamic contribution of helix residues in this domain. Together, these results provide a structural basis for understanding UBL3-dependent protein interactions and disease-related mechanisms, and suggest that helix-centered dynamic control in the UBL domain may represent a potential target for modulating UBL3 function.

Affinity chromatography is a powerful and therefore popular method for the purification of proteins for structural studies. The success of the technique relies on the specificity of the interaction between the target protein and the affinity resin. Here, we present the identification of two protein contaminants isolated from HEK293 cell lysate following affinity purification of twin Strep-tagged or FLAG-tagged proteins. The contaminants were identified as human propionyl-coenzyme A carboxylase (hPCC) and protein arginine methyltransferase 5 in complex with methylosome protein 50 (PRMT5:MEP50) via a combination of cryo-EM data processing and proteomic analyses. This report serves to illustrate how these contaminants may appear in cryo-EM datasets and to highlight the paramount importance of affinity chromatography resin specificity for efficient protein purification.

Our current understanding of carbohydrate-binding module (CBM) function is limited by the fact that most CBM research has focused on single-binding-site modules. CBM family 92 (CBM92) is a recently characterized family of predominantly trivalent proteins that bind {beta}-1,3- and {beta}-1,6-glucans with high specificity. CpCBM92A from Chitinophaga pinensis stands out as the first trivalent member of the family to be structurally determined. Multivalent CBM families are rare, and the way in which the three binding sites cooperate in ligand recognition remains unclear. Here, we use NMR spectroscopy to demonstrate how each of the proteins binding sites plays distinct roles in ligand binding. One binding site, referred to as the {beta} site, can be identified as the primary attachment point because of its higher affinity for all tested ligands, consistent with previous biochemical data suggesting it is the strongest binding site on CpCBM92A. The other two binding sites, referred to as and {gamma}, preferentially bind longer segments of {beta}-1,3- and {beta}-1,6-glucan chains, respectively. We further show that the glycosidic bond position and anomeric configuration of the binding glucosyl unit strongly affects protein affinity due to a preferred ligand pose in the binding sites. Our results provide insight into how the trivalent architecture of CBM92 might enable cross-linking of scleroglucan chains, which may guide the development of new applications for CBMs in biotechnology.

Human RNase 2 (eosinophil-derived neurotoxin, EDN) is a major eosinophil granule protein of the vertebrate-specific RNase A superfamily and is involved in antiviral response and inflammation. Identifying ligand-binding pockets in EDN is thus relevant to structure-based drug design. In our laboratory we identified by protein crystallography a conserved site at the protein surface binding to carboxylic anion molecules (malonate, tartrate and citrate). Searching for potential biomolecules rich in anion groups and considering previous report of EDN binding to glycosaminoglycans, we explored the protein binding to saccharides. Next, EDN crystals were soaked with mono- and disaccharides, and the 3D structures of ten complexes were solved by X-ray crystallography at atomic resolution. We identified protein binding pockets to glucose, fucose, mannose, sucrose, galactose, trehalose, N-acetyl-D-glucosamine, N-acetylmuramic acid, and the sialic acid N-acetylneuraminic acid. A main site for glucose, fucose, and galactose was located adjacent to the spotted carboxylic anion site. Secondarily, N-acetylneuraminic acid, N-acetylmuramic acid, sucrose, galactose, and mannose shared another protein surface region. Overall, the saccharides clustered into seven defined sites, outlining a conserved recognition pattern, which was further analysed by molecular modelling. Interestingly, within the RNase A family, we find amphibian RNases that were initially isolated as carbohydrate binding proteins and named as leczymes, combining enzymatic and lectin properties. The present data is the first systematic structural characterization of a mammalian sugar-binding RNase within the family. The results highlight unique EDN residues that mediate its sugar specific interactions, of particular interest for a better understanding of the protein physiological role. HighlightsO_LIstructure of RNase 2 in complex with mono and disaccharides at atomic resolution C_LIO_LIidentification of RNase 2 unique sugar binding sites C_LIO_LIcharacterization of a mammalian RNase A family enzyme with lectin properties C_LI Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=110 SRC="FIGDIR/small/713198v1_ufig1.gif" ALT="Figure 1"> View larger version (46K): org.highwire.dtl.DTLVardef@1d805f7org.highwire.dtl.DTLVardef@16fcc49org.highwire.dtl.DTLVardef@ccfd92org.highwire.dtl.DTLVardef@1b8f1e_HPS_FORMAT_FIGEXP M_FIG C_FIG

TBCK-related encephalopathy (TBCKE) is a neurodevelopmental disorder associated with biallelic mutations in TBCK. Despite the increasing number of reported cases worldwide, the biochemical and biophysical properties of TBCK remain unclear, hindering molecular understanding of its role in disease. Here, we present the successful expression, purification, and biochemical characterization of full-length human TBCK produced in Spodoptera frugiperda cells. Biochemical and biophysical analyses reveal that the catalytically inactive pseudokinase domain of TBCK lacks nucleotide binding, consistent with the absence of the canonical VAIK, HRD, and DFG motifs required for catalysis. These findings support that TBCK is a class I pseudokinase and provide a foundation for future structural and functional studies to elucidate its biological role.

FAM122A regulates cell cycle progression through inhibition of the PP2A-B55 phosphoprotein phosphatase. Recent structural work has uncovered helical elements in the N-terminus of FAM122A as binding determinants for PP2A-B55 but whether FAM122A inhibition towards PP2A-B55 is regulated is presently unclear. To address this we performed a systematic analysis of the PP2A-B55 interaction with FAM122A in cells uncovering a novel region in the C-terminus of FAM122A, spanning residues 150-170, required for binding. This C-terminal region and the N-terminal helices are both required for efficient binding to PP2A-B55 suggesting a bipartite binding mechanism. We perform amino acid resolution scans of FAM122A 150-170 uncovering several residues in this region contributing to binding including the conserved Ser158, a reported phosphorylation site. We show that Ser158 is important for PP2A-B55 inhibition in human cells as well as efficient stimulation of mitotic entry in Xenopus laevis egg extracts. In human cells and in Xenopus laevis Ser158 phosphorylation is regulated with increased occupancy correlating with cell cycle stages requiring PP2A-B55 inhibition. Collectively our work uncovers novel aspects of FAM122A interaction with PP2A-B55 and provides a possible mechanism for how the inhibitory activity of FAM122A can be regulated during the cell cycle.

Human immunodeficiency virus (HIV) infection promotes considerable bioenergetic, spatially heterogenous strain to the brain that is incompletely ameliorated through viral suppression afforded by antiretroviral therapy (ART). Disrupted homeostasis of brain lipids after HIV in humans or simian immunodeficiency virus (SIV) infection in rhesus macaques occurs due to elevated energetic demands, neuroinflammation, reactive oxygen species, and barrier leakiness. Brain lipids are particularly vulnerable to HIV-associated dysregulation due to their high abundance, unique composition, and specialized functional roles. Using rhesus macaques exposed to SIV and ART (tenofovir disoproxil fumarate (TDF), emtricitabine (FTC), and dolutegravir (DTG), we investigated the spatial distribution and abundance of lipids across brain regions and metabolically relevant peripheral tissues using mass spectrometry imaging. When comparing lipid abundance, individual lipids representing a multitude of species were more varied across tissues than by treatment condition. Further, we discerned either solely SIV infection or ART outweighed one another in altering phospholipids in different tissues Presence of ART had a greater influence on phospholipid homeostasis in the temporal cortex and hippocampus than in the midbrain, possibly due to differences in penetrance and turnover of ART across brain regions. Overall, these data demonstrate ART robustly increased phospholipids across brain regions while SIV infection had a varied impact depending on the brain region. These findings inform the need to further evaluate the neurologic consequences that may result in the brain due to disrupted lipid homeostasis across ART regimens.

Nicotinic acetylcholine receptors (nAChRs) belong to the pentameric ligand-gated ion channel superfamily (pLGICs). Among them, the neuronal homomeric 7 nAChR is highly permeable to calcium and plays critical roles in synaptic transmission, cell signaling, and inflammation modulation. The biogenesis of 7 nAChRs is enhanced by the chaperone proteins RIC-3 and NACHO. Previously, we reported a motif in the 5-HT3A receptor, another pLGIC, involved in RIC-3 modulation. Residues in this motif are conserved and also found within the L1-MX segment of the 7 nACh subunit. We therefore explored the regulatory roles of these conserved residues in the biogenesis of 7 nAChRs using multiple approaches, including heterologous expression in Xenopus laevis oocytes, mutagenesis, pull-down assays, cell-surface labeling, and two-electrode voltage-clamp (TEVC) recordings. We find that synthetic 7 L1-MX peptide interacts with both RIC-3 and NACHO. In particular, conserved residues W330, R332, and L336 in the L1-MX positively regulates the assembly of 7 oligomers and the biogenesis of 7nAChR. In presence of residues W330, R332, and L336, NACHO promotes an assembly of an 7 pentamer which is resistant to strong denaturing conditions. NACHO-promoted 7 pentamer is also resistant to Endo H enzyme. Sensitivity of the pentamer to moderate temperatures (37 {degrees}C, 45 {degrees}C, and 50 {degrees}C) suggests that NACHO stabilizes the pentamer via non-covalent interactions. In contrast, Ala replacements at these residues disrupt the biogenesis and abolish 7 current. NACHO and RIC-3 co-expression yields partial rescue of functional expression for some Ala replacement constructs. SUMMARYThis work identifies regulatory roles of conserved residues W330, R332, and L336 in the biogenesis of 7 nAChR. This discovery positions MX subdomain as a promising target for future drug development that can minimize adverse effects.

Glutamate racemase (MurI) catalyzes the stereochemical interconversion of L-glutamate to D-glutamate, a key element of bacterial peptidoglycan biosynthesis. In this study, we present the crystal structure of Helicobacter pylori glutamate racemase at 1.43 [A] and in monoclinic symmetry, as previously reported models, but different unit-cell parameters. The present model contains a single dimer and retains the previously described head-to-head dimer arrangement. The differences between the models arise from variations in unit-cell parameters, which lead to altered crystal packing interactions rather than changes in the quaternary assembly. The monomeric fold and active-site architecture remain conserved and are consistent with the catalytic features described for bacterial glutamate racemases. This structure provides an updated, high-resolution structural model for H. pylori glutamate racemase and highlights the variability in crystal packing within the same space group.

Calculation of density maps from atomic models is essential for structural studies using crystallography and electron cryo-microscopy (cryoEM). These maps serve various purposes, including atomic model building, refinement, visualization, and validation. However, accurately comparing model-calculated maps to experimental data poses challenges, particularly because the resolution of cryoEM experimental maps varies across the map. Traditional crystallography methods generate finite-resolution maps with uniform resolution throughout the unit cell volume, while most modern software in cryoEM employ Gaussian-like functions to generate these maps, which does not adequately account for atomic model parameters and resolution. Recent work by Urzhumtsev & Lunin (2022, IUCr Journal, 9, 728-734) introduces a novel method for computing atomic model maps that incorporate local resolution and can be expressed as analytically differentiable functions of all atomic parameters. This approach enhances the accuracy of matching atomic models to experimental maps. In this paper, we detail the implementation of this method in CCTBX and Phenix. SynopsisNew tools implemented in CCTBX and Phenix allow the calculation of variable-resolution maps through a sum of atomic images expressed as analytic functions of all atomic parameters, along with their associated local resolution.

GPI-anchored proteins are crucial cell surface proteins with diverse, organism-specific functions, in eukaryotes. They are produced when the GPI transamidase (GPIT), a five-subunit membrane-bound enzyme complex, attaches a pre-formed GPI anchor to the C-terminal end of nascent proteins on the lumenal face of the endoplasmic reticulum. This process requires the removal of a C-terminal signal sequence (SS) on the substrate protein by the action of an endopeptidase subunit of the GPIT, Gpi8/ PIG-K. Using an AMC-tagged peptide in a cell free (post-mitochondrial fraction) assay, this manuscript studies the steady state kinetics of enzymatic cleavage of the substrate by GPIT of the human pathogenic fungus, C. albicans. We show that Mn+2 enhances activity by improving substrate binding but plays no direct role in substrate cleavage per se. Molecular dynamics simulations suggest that the divalent cation binds at a site away from the active site but provides compactness and stability to Gpi8. It also enables a conformation in which a flexible loop (219-244 residues) in the vicinity of the catalytic pocket is able to interact with and position the scissile bond for cleavage by Cys202. Steady state kinetics also indicate that peptides of lengths 7-mer to 9-mer are better bound than 4-mer or 15-mer peptide substrates. A bulky residue at the site of cleavage reduces the catalytic activity of the GPIT. This is the first detailed steady state kinetics study on the endopeptidase activity of a GPIT from any organism.

There has been a surge in antibiotic resistance in recent years, making traditional antibiotics less effective against key pathogens. RNA has recently emerged as a potential target for antibiotics due to its involvement in crucial microbial functions. It is possible to expand the range of therapeutic targets by using RNA-based therapies, but it remains necessary to improve the molecular-level understanding of interactions between RNA and known and potential binders. The SAM-I riboswitch, which controls the transcriptional termination of gene expression involved in sulfur metabolism in most bacteria, is an excellent ligand target. Thus, understanding its behavior with and without ligand complexes would be very helpful for drug design applications. In this manuscript, we studied the interactions between the SAM-I riboswitch and its natural ligand, SAM, which controls riboswitch function, and compared those interactions to its interactions with the very similar small molecular SAH, which does not control riboswitch function, and to its interactions with a potential binder JS4, identified via virtual screening. From our simulations, we gain a deeper understanding of small molecule interactions with the SAM-I riboswitch. The results reveal how differently the small molecules (SAM, SAH and JS4) bind to and potentially induce conformational changes in the riboswitch. Our findings offer valuable insight into the molecular mechanisms underlying riboswitch RNA-ligand interactions for the design of more effective RNA-targeting therapeutics.

Toll-like receptors (TLRs) are vital components of the innate immune system, recognizing both exogenous pathogens signals (PAMPs) and internal stress signals (DAMPs). TLR2 is unique among the human (Homo sapiens) TLR family members, as it contains a large cavity for binding hydrophobic ligands, such as lipoteichoic acid (LTA) and di/triacyl lipopeptides (Pam2/3CSK4). This study analyzed the structural phylogeny of cavity presence in the TLR2 lineage in vertebrates (vTLR) enabled by AI protein structure predictions and explored the potential convergent evolution of similar features in invertebrates (iTLRs). Analysis of AI models of TLR2s shows that this cavity is consistently present in TRL2 orthologs across jawed vertebrates (Gnathostomata). In jawless vertebrates (Cyclostomatha), these cavities were found in lamprey (Petromyzon marinus) TLR2 model, but only in some extant hagfish (Myxini), suggesting an ancestral origin in basal vertebrates followed by lineage-specific losses. TLR2 paralogs were found in several species, with a similar central cavity but potentially different ligand specificities. In silico ligand docking showed Pam2CSK4 binds to this cavity in all TLRs and paralogs consistently, demonstrating the conserved function of the ligand-binding pocket in gram-positive bacteria recognition across TLR2 branches. Changes in the TLR2 cavity size and shape in some vertebrate groups show the evolution of this DAMP recognition mechanism adapted to its respective pathogens. iTLRs form a separate phylogenetic branch with distinct structural features, but in literature some are considered to be TLR2 orthologs. Indeed, TLRs from some species of Helobdella and Ciona, contain a cavity with some similarity to that in the vTLR2 lineage. However, detailed structural comparisons of their location in the LRR domain and the structural details of the models suggest that their cavities have developed independently from that in TLR2s. Smaller cavities are present in other branches of the LRR family, but show different locations, shapes, and features, indicating that the binding of small ligands in the internal cavities within the LRR domains evolved multiple times in the LRR domain family history.

Lacritin is an abundantly expressed glycoprotein in tear fluid and plays key roles in immune response, tear secretion, and bacterial killing. These biological functions are tightly regulated through several biochemical mechanisms including multimerization, proteolysis, and alternative splicing, especially within its C-terminal domain. Given its critical role at the ocular surface, lacritin is currently under investigation as a diagnostic biomarker and therapeutic candidate for dry eye disease (DED). However, despite over three decades since its initial discovery, the functional significance of the O-glycans that comprise more than 50% of its molecular weight remain largely unknown. To address this gap, we leveraged mass spectrometry (MS)-based glycoproteomics and molecular dynamics (MD) to explore the structural role of site-specific O-glycans on C-terminal lacritin. In doing do, we identified distinct glycosylation profiles between monomeric and multimeric lacritin, particularly at glycosites located near crosslinking residues (Lys101 and Lys104) that modulate multimer formation. Building on our glycoproteomics data, we performed MD simulations on monomer and multimer glycoforms and revealed that O-glycans participate in intra-glycan-protein interactions, thereby affecting the conformational flexibility of lacritin and the spatial arrangement of Lys101 and Lys104. Finally, we quantified the solvent-accessible surface area (SASA) of Lys101 and Lys104, highlighting that proximal O-glycosylation is predicted to affect the propensity of these residues to participate in crosslinking. Taken together, these findings underscore a central role for lacritin O-glycans in affecting structural topology with implications for its downstream biological activity.

Single molecule methods have become prevalent tools in elucidating molecular processes across various life science fields over the past three decades, driving breakthroughs in understanding their underlying molecular mechanisms. In our study, we employed two single-molecule methods, Forster Resonance Energy Transfer (smFRET) and high-resolution optical tweezers, to investigate the transcription of Mycobacterium tuberculosis RNA polymerase (MtbRNAP) from initiation through to termination. We aim to provide a set of comprehensive biophysical tools to deepen our current understanding of MtbRNAP and its transcription factors. These experimental assays represent an important step towards unraveling the molecular dynamics and interactions that support transcription in Mycobacterium tuberculosis.

Mechanosensing involves proteins detecting mechanical changes in the cytoskeleton or at cell adhesion sites. These interactions initiate signaling cascades that produce biochemical effects such as post-translational modifications or cytoskeletal rearrangements. Filamin is a ubiquitous mechanosensing protein that binds actin filaments and senses pulling forces within the cytoskeleton. Drosophila Filamin (Cheerio) is structurally similar to mammalian Filamin, with roles in egg chamber development, embryo cellularization, and integrity of muscle attachment sites and Z discs in Drosophila indirect flight muscles (IFMs). Here we report a potential novel binding partner of Drosophila Filamins: the death-associated protein kinase Drak that functions as a myosin light chain kinase. We found that Drak biochemically bound to an open mutant of Filamin that resembles the mechanically activated form partially bound to wild type Filamin and did not bind to closed mutant of Filamin. The interaction site was mapped to the intrinsically unfolded C-terminal region of Drak. To study the functional role of Drak-Filamin interaction, we studied two developmental events where Drak has been earlier shown to be expressed and where Filamin also functions: early embryonic cellularization and indirect flight muscle development at pupal stages. We found partial colocalization between Drak-GFP and Filamin-mCherry during the initiation of cellularization furrow, and at the time of myotube attachment site maturation in tendon cells. However, functionally we could not show direct correlation between Filamin and Drak. Our studies reveal interesting new expression patterns of Drak during Drosophila development and provide detailed information about Filamin localization during IFM development.

Synonymous codons encode the same amino acid but can differ in their usage and translational properties. In previous work we reported statistical differences in backbone dihedral angle distributions associated with synonymous codons in the Escherichia coli proteome. This finding has been questioned due to concerns regarding the statistical methodology used. Here we revisit the dataset using corrected statistical procedures and alternative statistical tests. Across multiple frameworks, the real dataset consistently shows an excess of small p-values relative to randomized controls, indicating detectable codon-associated differences in backbone conformation.

Protein translation is a highly regulated process influenced by multiple factors at the initiation, elongation, and termination stages. One notable regulatory element of the ribosome is the CAR interaction surface, a three-residue motif in the structure of the ribosome composed of C1274 and A1427 of S. cerevisiae 18S rRNA (corresponding to C1054 and A1196 in E. coli 16S rRNA) and R146 of ribosomal protein Rps3. CAR is highly conserved and positioned adjacent to the amino-acyl (A site) decoding center. It establishes hydrogen bonds with the +1 codon next in line to enter the ribosome A site, acting as an extension of the tRNA anticodon and forming base-stacking interactions with nucleotide 34 of the tRNA. However, despite CARs enzymatically strategic positioning within the ribosome, its functional relationship with the A site remains poorly characterized. Using molecular dynamics (MD) simulations, we examined the interplay between the A site and CAR site, revealing sequence-dependent modulation of H-bonding and {pi}-stacking interactions within and between the two sites. These findings highlight the interplay between the A site and CAR site, suggesting a structural and functional connection between these two regions of the ribosome that may contribute to mRNA sequence-specific tuning of translation elongation. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=91 SRC="FIGDIR/small/714784v1_ufig1.gif" ALT="Figure 1"> View larger version (22K): org.highwire.dtl.DTLVardef@1d783d3org.highwire.dtl.DTLVardef@f9cd8org.highwire.dtl.DTLVardef@102667corg.highwire.dtl.DTLVardef@967c56_HPS_FORMAT_FIGEXP M_FIG C_FIG

The cell surface localization of the Na,K-ATPase (sodium pump) is required for maintaining transmembrane electrochemical gradients. While glycosylation of the {beta}1 subunit facilitates trafficking from the endoplasmic reticulum to the plasma membrane, its role in nanoscale surface organization is not characterized. This study employed GlycoSHIELD computational modeling and DNA-PAINT single-molecule localization microscopy (SMLM) to evaluate how N-glycans influence pump distribution. In-silico simulations indicated that N-glycans sequester the protein core, providing a steric shield that increases with structural complexity. To investigate this experimentally, glycosylation-deficient mutants (3NQ) were generated and confirmed via immunoblotting. Quantitative SMLM analysis of A498 cells demonstrated that wild-type pumps exhibit higher localization density and form larger (144 nm) and more frequent clusters than 3NQ mutants (109 nm). These results indicate that N-glycosylation promotes stable enzyme clustering, supporting a galectin-lattice mechanism of organization rather than steric repulsion.

Synaptic vesicle glycoproteins 2 (SV2) are integral membrane proteins essential for neurotransmitter release and are implicated in neurological disorders including epilepsy and Parkinsons disease. In the attempt to develop a ligand selective for SV2C, and in collaboration with UCB, UCB-F was identified as a potential candidate. However, the affinity of UCB-F to SV2C was found to be temperature dependent, decreasing by about 10-fold from +4 to 37 degrees. UCB1A was subsequently identified as SV2C ligand displaying in vitro a 100-fold selectivity for SV2C compared with SV2A. In this study we investigated whether the binding of UCB-1A to SV2A and SV2C was affected by the temperature. A combination of experimental binding assay data and molecular dynamics (MD) simulations were used. The binding studies revealed that UCB1A affinity for SV2A decreased significantly at 37 {degrees}C compared with 4 {degrees}C, whereas binding to SV2C remained largely unchanged. MD simulations reproduced these observations, namely that ligand RMSD values at 310 K showed that UCB1A binding fluctuated markedly in the SV2A complex, with many trajectories exceeding the 3.0 [A] stability cutoff, whereas UCB1A remained relatively well-anchored in SV2C under the same conditions. Structural analysis showed that, while UCB1A adopts a conserved binding pose across all isoforms stabilized by {pi}- {pi} stacking and a hydrogen bond with Asp, SV2C possesses a unique stabilizing feature. In SV2C, Tyr298 is less exposed to the solvent and engages in a persistent hydrogen bond with Asparagine, a structural feature that reinforces pocket stability and limits temperature-induced destabilization. This interaction is absent in SV2A, consistent with its greater temperature sensitivity. Together, these findings provide a mechanistic explanation for the experimentally observed temperature independence of UCB1A binding to SV2C. More broadly, the results highlight the importance of incorporating physiologically relevant temperatures into SV2 ligand evaluation and demonstrate how combining experiments with simulations can uncover isoform-specific mechanisms of ligand recognition and stability.